Sensitive Chemical Compass Assisted by Quantum Criticality 
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The radical-pair-based chemical reaction could be used by birds for the navigation via the geo- 
magnetic direction. An inherent physical mechanism is that the quantum coherent transition from 
a singlet state to triplet states of the radical pair could response to the weak magnetic field and 
be sensitive to the direction of such a field and then results in different photopigments in the avian 
eyes to be sensed. Here, we propose a quantum bionic setup for the ultra-sensitive probe of a weak 
magnetic field based on the quantum phase transition of the environments of the two electrons in 
the radical pair. We prove that the yield of the chemical products via the recombination from the 
singlet state is determined by the Loschmidt echo of the environments with interacting nuclear spins. 
Thus quantum criticality of environments could enhance the sensitivity of the detection of the weak 
magnetic field. 
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Introduction. — Since Schrodinger questioned "what is 
life" from the general point view of a quantum physi- 
cist [l[ , scientists have never stopped the long-term explo- 
ration for the physical sources of the living phenomena, 
and this even stimulated the enthusiasm for the great 
discovery of the DNA genetic molecule [2^. Today it 
seems trivial to say that the life is of quantum for the 
molecules composing lives obey quantum laws, but some 
recent discoveries are very intriguing as some optimized 
living processes may be based on the nontrivial quantum 
effect from quantum coherence. One example in point 
is the photosynthesis process. Recent experiments have 
been able to exactly determine the time scales of various 
transfer processes by the 2D optical spectroscopy, and 
then show quantum coherence effects in energy transfer 
via collective excitations of some light-harvesting com- 
plexes [3|-[s|- 

Another prototype of quantum coherence effect for liv- 
ing process seems to appear in the avian magnetorecep- 
tion mechanisms [gI, verified by some recent experi- 
ments [8]. Recently quantum information approaches 
have been used to further analyze the role of quantum 
coherence phenomenon in the avian magnetoreception 
models [9-12]. It is now believed [Ts*, Ti] that the model 
for magnetoreception is based on the radical-pair mech- 
anism (RPM) [6] : the radical-pair molecule with two un- 
paired electrons is activated by light. When the electrons 
interact with their individual environments of nuclei via 
the hyperfine couplings, the spin singlet state will transit 
to the spin triplet states even though the external field 
is uniform and rather weak. In response to this quan- 
tum coherent transition, the field-dependent change in 
the product yield of the radical-pair-based chemical re- 
action is enough to be sensed by the avian retina. 

In spite of the rapid progress in the understanding of 
the RPM in the last decade, a key question remains 
elusive [13]: why does the singlet-triplet interconver- 




Radical Pair 

FIG. 1: (color online). Model setup of sensitive magneto- 
detection based on the radical pair chemical reaction assisted 
by the quantum critical environments consisting of interact- 
ing spins. Two spins could be initially prepared in the sin- 
glet state, and finally evolve into the triplet states due to the 
couplings to their environments. The corresponding chemical 
product sensitively responses to the external magnetic field. 



sion response to the extremely-weak geomagnetic field 
10~^T), and why is it very sensitive to its direc- 
tion? We notice that the existence of nuclear environ- 
ments surrounding the electron spins in the radical-pair 
molecule is crucial to magnetic sensitivity of the chemi- 
cal reaction. This observation motivates us to consider 
the role of internal quantum correlation in each envi- 
ronment. In this Letter, we propose a quantum phase 
transition (QPT)-assisted setup for the probe of a weak 
magnetic field in understanding the above conundrum. 
Actually a lot of real-world detectors are built based on 
the dramatic changes of the systems around phase tran- 
sitions, which amplifies the ultra- weak signal and thus 
enable one to probe it. Examples include bubble cham- 
ber detectors [l5| and superconducting single-photon de- 
tectors [16], where the liquid-gas phase transition and 
superconductor-metal phase transition take place respec- 
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lively enhancing the sensitivity for detection. We calcu- 
late the corresponding chemical product yield, which is 
phenomenologically described by a damping process [l^ • 
It is discovered that the chemical product yield is deter- 
mined by the time integral of the Loschmidt echo (LE). 
Our result shows that the chemical compass assisted by 
quantum criticality indeed can response to a very-weak 
magnetic field and be sensitive to its direction. 

Quantum- Phase- Transition- Assisted Radical Pair 
Mechanism. — Our setup is illustrated in Fig.[Tl Each of 
the two electrons in the radical pair is uniformly coupled 
to its own quantum correlated environment, which can be 
described by a transverse field Ising (TFI) model isl, [HI . 
In an external magnetic field B {cos Ox + sin^^), which 
may be the geomagnetic field, the environment is 
described by H'^ = + qn/^nB sinOI^ j with 
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Hn — J / ^ (-^n,.7^n,.7 



/ ^V"-n,.7--n,.7 + l 



(1) 



where A = qniinB cosO / J is the rescaled strength of 
the transverse field in unit of J being the Ising coupling 
constant, QnIJ^n is the nuclear magnetic moment, and 
n = 1,2 refers to the environment of the nth electron. 
I^ j and I^ j are Pauh matrices of the jth nuclear spin 
operators. In case of antiferromagnetic Ising chain, i.e., 
J > 0, we can omit the longitudinal terms and H'^ ^ Hn 
since they lead to only higher-order correction |2fl . Here, 
different from all previous studies 0, EH, 0,0, we 
explicitly consider the inter-nucleus couplings ^n,j^n,j+i- 
This coupling competes with the Zeeman energy being 
proportional to the geomagnetic field, and leads to a QPT 
at the critical point A = 1. A central spin uniformly 
coupled with each spin in this TFI system possesses the 
dynamic sensitivity described by the sharp decay of LE 
near the QPT [22|, which has been experimentally veri- 
fied [13, 

The two unpaired electrons in the radical pair couple 
to the two environments Ei and E2 respectively with the 
following Hamiltonians 



Vn = ^ sin l9cr^ + Vt cos 6>cr^ + Jgcrl^ 



(2) 



where and for n = 1 , 2 are the Pauli operators for 
the nth electron spin, the dimensionless coupling con- 
stant scales as ^ = gQ/\fN in the van Hove limit for 
the interacting many-body system. All the information 
about the geomagnetic field is also incorporated in and 
r^, which is the electronic Zeeman energy splitting in- 
duced by the geomagnetic field. 

Radical-Pair Evolution in Correlated Nuclear 
Environment. — Due to the spin-fiip terms, the 
time evolution governed by the total Hamiltonian 
H = ^n^Hn + Vn) can only be solved with some approx- 
imation. Usually the electron spins evolve faster than the 



nuclear spins. Thus we can first regard nuclear spins as c- 
numbers for formally diagonalizing the electronic Hamil- 
tonian Vn through a generalized Born-Oppenheimer 
approximation [25]. The eigen states of the electron 
spins are obtained as |+) = cos(a/2) W) + sin(a/2) |^) 
and |-) = sin(a/2)|t) - cos(a/2) ||) where |t) 
(1^)) is the spin- up (down) state in the Gx- 
representation with corresponding eigen val- 
ues E± = ±E for E = ^Q? sin^ l9 + and 
A = 1] cos -\- Jg . And the mixing angle is 

defined as a = 7r/2 — tan~-^ {ftsmO/A) . 

For the weak coupling {J go <C 0.) of an electron to 
nuclei we shall approximately obtain the eigen states 
1+) c^cos6>'|t) + sinl9'|;) and |-) sin |t) - cos |;) 
to the zeroth order of g for 0^ = {7r/2 — 6) /2 and the 
eigen energy E ft -\- Jg cos 6 to the first order. 

The Born-Oppenheimer approximation shows that the 
slowly- varying nuclear spins would not induce the coher- 
ent transition of the fast- varying electronic degrees, but 
the electronic motion provides an effective potential for 
the nuclear spins. In this sense, the total Hamiltonian is 
approximately rewritten as 



(3) 



where the different effective Hamiltonians for the nuclear 
spins corresponding to states |±) are respectively 

^" = J2j Vnjn,j+1 + (A ± 5 COS 0)11^^] ± n. (4) 

It can be proven that the Born-Oppenheimer approxima- 
tion is generally valid even for such a large N that a QPT 
occurs. 

Product Yield and Loschmidt Echo. — The radical pair 
is assumed to be initially in the singlet state l^*) which 
subsequently suffers from the homogeneous interaction 
V = Vn with the environmental nuclear spins. Then 
the radical pair undergoes a singlet-to-triplet transi- 
tion. The charge recombination of the radical pair goes 
through different channels, depending on the electron- 
spin state (singlet or triplet). In particular, the singlet- 
state product yield formed by the reaction of radical pairs 
can be calculated as [l7| 



f 

Jo 



rcit)fsit)dt, 



(5) 



where rc{t) is the radical re-encounter probability distri- 
bution, and fs{t) = {S \pe{t) \ S) is the singlet-state pop- 
ulation at time t. Usually it is assumed that rdt) = 
ks {—kst) with ks the recombination rate. The ulti- 
mate product yield = ^s{t 00) in cryptochrome 
is believed to affect the visual function of animals [6] . In 
order to quantitatively describe the magnetic-field sen- 
sitivity of the radical-pair reaction, we shall resort to 
A(^) = 9$5'/96>,the derivative of the product yield with 
respect to the geomagnetic-field direction [2L] . 
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For nuclear spins initially in the mixed state pi (8) P2, 
the initial state of the total system is p(0) = liS) (iSI (g) 
Pi P2- When the state of the electron is |+) or | — ), 
the evolution of the environment is governed by or 
H~ respectively, and hence the initial state pn will evolve 
into (^ri^)^ or U~pn{U~)^ respectively. Here, the 

evolution operators are defined as = ex.p {—iH^t). 
This will result in the so-called adiabatic entanglement 
between the system and the environment [25]. For the 
above initial state p(0), at time t, the reduced density 
matrix for the electron spins pe{t) =tTn ('^(^)l reads 



To the opposite, we consider the case with a small A^, 
in which the above analysis fails. Therefore, we shall deal 
with it separately. For A/" = 2, we obtain explicitly 



L{t) = 1 



1652cos2 6'sin2(Ct) 
[l + 4(A-5cos6')2]^2 



(11) 



where ^ = ^/l -\- 4(A + gcosO)'^. It shows that the setup 
for a small N can not work as well as it for a large N. 
The detailed discussions about dependence of on N 
are visually given in the supplementary material. 



Peit) = -[| + ->{ + -| + |- + )(- + | 

-D{t)\+-){-+\-D*{t)\-+){+-\], (6) 

where D{t) = tr [Z+pi (C/f )^ tr p2 (U:^)^ is the 

P2 = p, D{t) is real and 



decoherence factor. When pi 



tr 



which is 



can be simplified to L{t) 
just the Loschmidt Echo characterizing the dynamic sen- 
sitivity of the environment in response to the perturba- 
tion [22]. Straightforwardly, we can prove the population 
fs{t) = [1 + /2, and obtain the product yield as 



1 1 



L{t)e- 



-kst 



dt. 



(7) 



This is a central result of our paper which reveals the 
direct relationship between the product yield and the LE. 

Hereafter, we calculate the product yield for identi- 
cal environments in an initial pure state, which cor- 
responds to the case with absolute zero temperature. 
In this case, the initial state can be described by 
a state vector \G) and LE is simplified as L{t) = 
\{G\exp{iH-t)exp{-iHn) |G)| , which was explicitly 
given in Ref. [22]. In the following we take two specific 
cases into account. First, we consider the large- case 
since the QPT usually occurs in this limit. Using the 



analytic results about L{t) obtained in Ref. [22|, we ap- 
proximate the product yield around the critical point as 




erf 



ks Y 
2^/7/ 



where erf (x) is the error function, and 
SJ^g^NK^ cos2 
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37r(l - A)2 



(8) 



(9) 



with Kc the cutoff momentum. For a sufficiently small 
i.e., ks <C 2y^, the product yield is approximated as 
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1 7rks\l-X\ 

2 WJgcosO V NK^' 



(10) 



Since d^s/dX is discontinuous when A = 1, it may serve 
as a witness of the QPT. 




FIG. 2: (color online). The product yield $s vs the mag- 
nitude B and direction of the magnetic field at a finite 
temperature T — 0.2 with N — 1000, go — 1, ks — 0.1, and 
J= 1. 

Sensitive Magnet odetection at Finite Temperatures. — 
The above results are obtained for an ideal case with pure 
states. For the practical purpose, we need to consider the 
cases at finite temperatures. In order to illustrate the re- 
sult for a very-large A/", we numerically plot the product 
yield <I>5' vs the magnitude B and direction of magnetic 
field at a finite temperature in Fig. [2l Obviously, the 
product yield displays its dependence on both the geo- 
magnetic field's magnitude and direction. Besides, there 
is a deep valley around the top- left corner. This can be 
seen from the fact that the LE decays in a gaussian way 
around the critical point A ^ 1. Additionally, in the re- 
gions far away from the critical point, e.g., at the top- left 
and bottom-right corners, the product yield nearly stays 
unity for the LE almost does not decay. 

Furthermore, in order to investigate the infiuence of 
other parameters, we plot the product yield vs direc- 
tion for different temperatures and recombination rates 
in Fig. [3l The similarity among the cases with different 
temperatures is that there is a peak for A(^) as it in- 
creases from zero at ^ = 0. It is seen that as the temper- 
ature increases, the position of the peak moves towards 
= 7r/2, meanwhile the line shape on the left hand side 
becomes more and more flat. In the high-temperature 
limit, we would expect a sharp peak around = 7r/2, 
while there is a platform elsewhere. In this case, the bird 
can no longer discriminate the direction. This is a rea- 
sonable result since a QPT takes place at the absolute 
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zero, and a high temperature sraears the QPT. When 
we come to the recombination rate ks in Fig. [3l^c), we 
see that the slower the singlet state reacts, the more the 
product yield changes along with the direction. That is 
because a reaction with a smaller recombination rate pro- 
vides more time for the decay of LE. The magnetic-field 
sensitivity in Fig. [H^d) clearly confirms our analysis. In 
addition, as the environment involves more nuclear spins, 
the visibility rises as the LE decays faster for a larger N. 



to other models as well, such as the so-called long-range 
TFI model, or Lipkin-Meshkov- Click model [19]. De- 
tailed studies of these generalization will be given in a 
forthcoming paper. 
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FIG. 3: (color online). The product yield and its deriva- 
tive A vs the geomagnetic field's direction for different T 
in (a) and (b) with ks = 0.1, and for different ks in (c) and 
(d) with T = 0.01. In all figures, we set N = 1000, go = 1, 
B 0.9 and J = 1. 

Conclusion. — We proposed a RPM-based magnetode- 
tection scheme assisted by the QPT. We have proved that 
the yield of the chemical product via the recombination 
from the singlet state is determined by the LE of the en- 
vironments. This relation results in the enhanced sensi- 
tivity of the RPM-based avian compass. Thus, our study 
does not only provide important insights to the mecha- 
nism of magnetoreceptor through a radical-pair process 
in a very- weak field, but it also shields new light on build- 
ing a quantum bionic device for ultra-sensitive magnetic- 
field sensing. In addition, in our bionic setup, the sen- 
sitivity is pronounced when the nuclear spin number is 
large and the recombination rate is small. 

It may be argued that in the avian retina, the envi- 
ronments of radical pairs consist of a few nucleus rather 
than numerous nucleus, i.e.. A/" ^ oo for a QPT. But 
the experiments 20|, 23, 24| however demonstrated that 
even for A^ = 2, there still exists the dynamic sensitiv- 
ity induced by quantum criticality. This result implies 
that dynamic sensitivity may have a close relation with 
the level crossing. Besides, although our scheme requires 
a very-low temperature for the current experimental pa- 
rameters, its importance also lies in the possible bionic 
setup for sensitive magnetodetection. Last but not least, 
our results are based on the TFI model, but the en- 
hancement of LE decay due to QPT is independent of 
the model [26]. Thus it is reasonable to infer that all 
the results obtained from TFI model can be generalized 
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